Accepted to ApJL 

Preprint typeset using I^TgX style emulateapj v. 6/22/04 



O 

o 



CONTINUUM ACCELERATION OF BLACK HOLE WINDS 

John E. Everett and David R. Ballantyne 

Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada 

Accepted to ApJL 

ABSTRACT 

Motivated by recent observations of high-velocity, highly ionized winds in several QSOs, models of purely 
continuum-driven winds launched from ^ 200 GA/bh /(? are presented. Launching conditions are investi- 
gated, as well as the observational signatures for a variety of initial conditions and illuminating continua. While 
we verify that continuum-driven, highly-ionized outflows reach the observed velocities for L/L^dd > 1, in- 
dependent of the incident spectral shape, such winds are too highly ionized to exhibit the observed absorption 
features when launched with an AGN continuum (in fact, such winds are so ionized that they are driven pri- 
marily by electron scattering). If the wind is instead illuminated with a blackbody continuum originating from 
an optically thick shield, the gas is too weakly ionized and does not produce high energy absorption features. 
If high-velocity, high-ionization winds are truly launched from very near the black hole, such winds must be 
launched under other conditions or via other processes; we summarize some possibilities. 
Subject headings: galaxies:active - quasarsiabsorption lines - accretion, accretion disks - galaxies:individual 
(PG 121lH-143)-X-rays:galaxies 
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1. INTRODUCTION 

Optical, UV, and X-ray absorption troughs in Active Galac- 
tic Nuclei (AGN) spectra are central to the ongoing ef- 
fort to understand not only mass outflow in AGN, but to 
provide clues to the fueling of their central superm assive 
black holes (see ICrenshaw. Kraemer. & Georgell200l . Al- 
though the optical depth, metallicity, and velocity of such 
outflows are becoming increasingly well-constrained, the 
mechanism accelerating the winds remains elusive: thermal 
(e.g.,EK i'olik & Kriss 1995), continuum-driven (for instance, 
Oieiouche & Netzer 2003), line-driven (e.g.. iMurrav et af] 
[1995; Proga. Stone. & Kallman 2000) , and hyd romagnetic 
(e.g., Konigl & Kaitje 1994; Bottorff et al]ll997h winds may 
all provide insight, but none of the above models have 
been excluded or are definitively observationally preferred for 
AGN outflows. 

In this Letter, observational constraints on continuum- 
driven winds launched from very near the central black 
hole are examined. Such wind models have become sig- 
nificant, as observations of outflows apparently moving 
at large fractions of the spe ed o f light (v ~ 0.1- 
0.4c) h ave been publ ished bv ICh artas. Brandt & Gallagher 
(|2003), Chartas et al.' ('2003'), Pounds et al. (2003a b), 'and 
iRee ve s, O'B rien. & Ward (2003). From the measured high 
velocities, the winds a re inferred to be launch ed very close to 
the central black hole. iKing & Poundsl (120031:') have proposed 
that the highly-ionized gas is subject to continuum driving 
near the central black hole because line driving is unlikely 
to be important for the large ionization parameters in these 
winds. 

A model of a continuum-driven wind is presented in S|2] 
This model is used to investigate terminal velocities and ab- 
sorption features displayed by a variety of winds with differ- 
ent incident continua, initial shielding columns, and densities 
at the base of the wind (fj3}- 
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Electronic address: [everett, ballantyne@cita.utoronto.ca| 



Our simple continuum-driven wind model consists of a ra- 
dial streamline set to rise above the disk at a constant angle 
of 5° to the disk, with the initial column (iVH,o) and num- 
ber density (uq) of the wind specified by the user. Gas with 
that initial density is assumed to be in Keplerian motion and 
is perturbed into the central, illuminating continuum with an 
initial velocity of roughly the local sound speed. To deter- 
mine the gas's motion due to continuum acceleration. Cloudy 
(Ferland 2002) simulations are run on a radial rays through 
the wind, spaced logarithmically in co-latitude. The contin- 
uum acceleration on the gas is then computed using both the 
ionization state and transmitted continuum (therefore taking 
account of any self-shielding effects) at the end of each radial 
slice through the wind, taking into account both electron scat- 
tering and bound-free opacities. The radiative acceleration is 
calculated and then compared to gravity in Fcont, defined as: 

p _ acont _ j-cIxM,^ 

where Ocont is the continuum radiative acceleration, g is the 
gravitational acceleration, p is the density in the wind, c is 
the speed of light, Xf is the frequency-dependent continuum 
opacity (output by Cloudy), is the transmitted continuum 
(again, given by Cloudy), G is the gravitational constant, 
i\/BH is the mass of the central black hole, and R is the spher- 
ical radius. This Fcont factor is then input into the Euler equa- 
tion along with the component of gravity and centrifugal ac- 
celeration along the radial flowline: 

dv GAfBH(l- Fcont) 



dr 



i?2 



Of) 



sin 6t, 



(2) 

where m is the cylindrical radius, 9 is the co-latitude, 9p 
is the (fixed) angle of the flowline from the vertical, and 
Vci,{m) — florQ/ru is the azimuthal velocity (angular momen- 
tum is conserved in outflowing gas). 

We integrate this equation to compute the velocity of the 
wind as a function of distance along the radial streamline. 
As that velocity increases, the number density in the wind 
(as well as the column density) decreases via mass conserva- 
tion. For the first set of Cloudy simulations, where no velocity 
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data are available, the density is set to fall as r^^; this initial 
density distribution is replaced in subsequent iterations and is 
unimportant for the final results. The solution is iterated at 
least twice to ensure that the density structure is consistent 
with the velocity profile in the wind; simulations in test cases 
have shown that the wind structure does not change after two 
iterations, even when ten iterations are completed. 

The wind parameters chosen in this paper reproduce 
the observational c onstraints on the outflow observed by 
iPounds et alJ (|2003a) in the narrow emission line quasar 
PG 1211+143. They observed an outflow with v ^ 0.08c 
and an inferred density of n ^ 3 x 10* cm~'^, situated 
at 3 X 10^^ cm from the black hole (this density is esti- 
mated observationally from t he emission measu re and line- 
of-sight column density; see Pounds et al Mm^ . They also 
inferred a launching radius of tq = 10^^ cm by noting that 
V = 0.08c ^ Vcsc, the escape velocity for that tq. Since 
M = Airhr^ pv (where h is the covering factor of the flow, r 
is the radial distance, p is the mass density, and v the veloc- 
ity of the wind) must be constant along the flow, an original 
density of 5 x lO^'^ cm^^ is inferred for the base of the 
wind if the initial vel ocity is of order 20 k m/s and b — 0.8 
(the value of b used in lPounds et alJ 2003aV Observations in- 
dicate that Mbh ~ 4 x 10^ Mq (Kapsietal. 2000)', and 
^Edd w 5 X 10"*^ erg/s (L/LEdd ~ 0.8 as the "Big Blue 
Bump" that dominates the continuum in PG 1211+143 has 
L ^ 4 X IC*^ erg/s; Pounds et al. 2003a|). The infeiTed mass 
outflow rate from the observations is ^ 2.8 Mq/ji (for the 
covering fraction b = 0.8). We assume solar abundances. 

The increase of density towards the base of the wind (where 
the gas is outflowing much more slowly) results in a large 
increase in the column density over the observed column o f 
~ 5 X 10^^ cm^^. As recognized bv lKing & Poundsl j2003h . 
this column would be significantly optically thick, preclud- 
ing photoionization simulations. Two possibilities are there- 
fore considered in our simulations. First, "standard" winds 
are tested with A'h.o < 10^^ cm~^ at the base of the outflow 
and with an AGN spectral energy distribution (SED). Next, 
winds are launched by a soft X-ray/UV blackbody from the 
photosphere of an optically thick shield, as hypothesized by 
lKing& Pounds (2003). 

3. KINEMATICS AND ABSORPTION FEATURES 

The most immediate check for these outflows is to see 
which parameters lead to the observed velocities. Figure 
displays Voo as a function of the Eddington ratio (voc here is 
the velocity at the end of the model's streamline at a distance 
of lOOro where rg is the launch radius of our wind; taking 
the terminal velocity at 10^ Tq only changes v^c by ~ 1% in 
tests). This figure shows the variation in Voo for both an AGN 
SED (Risaliti & Elvis 2004)- and for a blackbody continuum 
with a range of temperatures (as proposed by King & Poundsl 
120031) . As a simple starting case, the winds presented here are 
optically thin, with A^h.o — 10^^ cm^^; higher columns are 
considered below. 

As expected, when L/ L-^dd < 1, the wind does not reach 
^'obs ^ Vcsc, even with the added bound-free opacity used to 
calculating the continuum acceleration. For L/LEdd > 1, 

' A recent reanalysis bv 'Peterso n et alJ 12004) yields 1.46 ± 0.44 X 
10* Mq (with significant uncertainty noted in the analysis); we retain the 
value of 4 X 10^ Mq to compare with previous models. 

- The Risaliti & Elvis 1 2004) spectrum is approximated by using Cloudy 's 
"AGN" continuum with the Big Blue Bump peaking at T = 1.5 X 10^ K, 
OQX = -1-43, auv = —0.44, ctx = —0.9. 



winds can achieve the observed velocities. Figure ^ clearly 
shows that the blackbody continuum can also accelerate a 
wind, r esulting in almost identical velocities as when the 
iRisaliti & Elvis (2004) SED is used. This seems reasonable 
since the majority of Fcont is given by frequency-independent 
electron scattering. 

Next, to test the sensitivity of these results to initial param- 
eters, both the column and the initial density of the wind are 
varie d (see F ig. 13 for winds illuminated with the AGN SED 
from' Risahti &Elvi3 (12004 . This shows that changing the 
initial density by an order of magnitude (which increases M 
by the same factor) does not change the velocity structure of 
the wind significantly, due to the fact that the continuum is 
strong enough to highly ionize both winds. For such highly- 
ionized continuum driven winds, the degree of acceleration is 
governed by / p which is nearly identical for both winds (al- 
though the higher-density wind does have a lower ionization 
parameter near the base of the wind, yielding larger accelera- 
tions there than in the lower-density wind). 

Also in Figure 12] the shielding column at the base of the 
wind is varied from the normal A^h.o = 10^^ cm^^ to 
^H,o = 5 X lO^'^ cm^^ to investigate the "self-shielding" 
of a thicker wind. This thicker column also does not change 
the final velocities in the outflow. 

d I — • • • • 1 • • ' • ' ' • • ' 1 • • ' • 1 

- Risaliti & Elvis (2004) AGN SED 

- — — — Blackbody, T=10^ K 

- Blackbody, T=10*-^ K 
^ Blackbody, T=10* K 
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Fig. 1 . — The final velocity of continuum-diiven wind models as a function 
of the Eddington ratio of the central continuum for a range of incident SEDs. 
For continuum driving in such highly ionized sources, the final velocity is 
relatively independent of the SED. The observed velocity of v^hs = 0.08c 
in PG 1211+143 is shown for comparison. Note that, for i/-LEdd ^ !> 
winds do not escape the system, but we plot their velocities at s = lOOro to 
show the trend as L/I/Edd increases. 



These figures address the range of Eddington ratios neces- 
sary to launch a wind to the observed velocity. But what are 
the particular spectral signatures of such a wind? To answer 
this question. Figures |3] and 0] show the X-ray ionization pa- 
rameter, ^, as a function of L /L-^dd for those different mod- 
els. We define f in the same manner as Cloudy and XSTAR: 
47ri^x/?iH, where Fx is the integrated X-ray flux from 1 to 
1000 Rydbergs. To pick one 'fiducial' ^ out of the range of 
ionization parameters in the accelerating wind, ^ is extracted 
at a distance along the flow of s = 3ro (chosen to compare 
with the inferred distance of the wind in PG 121 1+143). 

Figures Inland 13 display the ionization parameters in these 
continuum-driven winds. They show that the high velocity, 
continuum-driven winds must by necessity have been accel- 
erated by a large incident flux, which, for AGN-like continua, 
leads to extremely high ionization parameters in the wind (^ 
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Fig. 2. — The final ve l ocity o f continuum-driven wind models vs. L/L^^^ 
for the lRisaliti & Elvil 120041) incident spectrum with changes to the initial 
density and initial wind column. Unless otherwise noted in the legend, the 
models have A^h.o = lO^'^ cm~^ and no = 5 X 10^^ cm~'^. Neither the 
increased initial column nor density affect the wind significantly. The ob- 
served velocity of fobs = 0.08c in PG 1211-H143 is shown for comparison. 



■Risaliti & Elvis AGN SED 

Risaliti & Elvis AGN SED, n„ = 5xlo" cm"^ - 
Risaliti & Elvis AGN SED, N„ = 5x10^^ cm"^ - 
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Fig. 4. — As in Fig. ^ but using the Risaliti & Elvis i2004) contin- 
uum to investigate variations in density and shielding. Unless otherwise 
noted in the legend, the models have A^jj o = 10 cm~^ and no = 
5 X 10^^ cm~'^. The larger density wind yields a somewhat smaller ^, even 
though it experiences a very similar acceleration to the lower-density model. 
The A'^H.o = 5 X 10^^ cm~^ wind has a very similar ionization parameter 
for L/iEdd > 1- 



of order 10^ - 10^ for i/Ledd > !)■ The T = 10^ K black- 
body spectrum also has a very high ionization parameter, but 
all of its energy in the 1-1000 Rydberg window is concen- 
trated near 1 Rydberg. Thus, even though this continuum has 
a large ionization parameter, the lack of X-ray photons means 
that the ionization stages within the gas are much lower than 
that produced by an AGN SED: for instance, the predomi- 
nant state of iron for the wind illuminated by the T — 10'^ K 
blackbody spectrum is Fe IX, while the obs erved iron ions are 
Fe XXV and Fe XXVI JPounds et alJ2003ah . 



Risaliti & Elvis (2004) AGN SED 

— — — Blackbody, T=10^ K 

Blackbody, T=10''-^ K 

Blackbody, T=10^ K 
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Fig. 3. — The X-ray ionization parameter, as a function of L/L-^dd 
for a range of incident continua at the inferred location of the absorbing gas 
in PG 1211. The relative values of the ionization parameter are set by the 
amount of flux each SED has in the range of 1 to 1000 Rydbergs, where § is 
calculated. The AGN SED yields ionization parameters too high to explain 
the observed line absorption. While the blackbody continua can have rela- 
tively high ionization parameters, they have no flux in the X-ray regime to 
ionize iron to the states observed in PG 1211-1-143. 



This has important implications for observations of such 
winds. King & Pound^ i2003h postulated that contin- 
uum driven winds could account for high-velocity X- 
ray absorption features. As mentioned in the in- 
troduction, there have been other recent observations 
of such high-velocity X-ray absorp tion i n a variety 
of QSOs jChartas. Brandt & Gallagher l2003t IChartas et all 



12001 IPounds et alJl2003bt IReeves. O'Brien. & Wardll2()0l . 

In observations of PG 1211-1-143, significant equivalent 
widths (EWs) for absorption features in Mg XII, S XVI, 
Fe XXV, and Fe XXVI are seen. Table^compares these obser- 
vations to models of a high-column wind with an AGN SED 
jRisaUti & Elvisl2004l l). iVn.o = 5 x lO^^ cm , assumed line 
widths of 1000 km/s ([Pounds et al. 2003^ and L/LEdd = 1- 
Equivalent widths are calculated in the linear weak line limit: 

Wcv = ASD\/^roX [l-ro/(2!%/2)+ToV(3!V3) + ...] (3) 

for ''0 ^ 1. where tq is the optical depth in the line from 
Cloudy and AE'd is the Doppler width of the line in eV (all 
of the lines shown here have tq < 1 in the wind). A com- 
parison of observed EWs to the model EWs clearly shows 
that, because of the high ionization in the outflow, such winds 
cannot account for the EWs observed in PG 1211-1-143. The 
EWs for winds illuminated with a blackbody spectrum are 
even smaller, although in that case the small EWs are due to 
the lack of high-energy ionizing X-rays. 

To te st the particular parameters presented bv lPounds et alJ 
(|2p03a), we assembled a model with a constant column of 
JVh = 5 X 10^"^ cm^^ present throughout the wind. Although 
the density still declines with height above the disk as before, 
we set the column to a constant value for every line of sight 
through the wind. This would therefore yield t he column of 
iVn = 5 X IQ-^'^ cm~^ at s = 3ro as observed bv'Pound s et alJ 
(2003a). This model differs strongly from the other solutions 
where the column drops consistently as the density in the wind 
falls. The last column of Tabled displays the EWs from this 
model with L/L-^dd = 1; at higher L/L-^dd, the EWs only 
decrease. Even with the constant column supplying a greater 
column of the relevant ions, the ionization parameter in the 
model is so high that the line equivalent widths fall below the 
90% confidence interval of the measurements. 

We have searched for some area of parameter space that 
might allow such continuum-driven winds to show the ob- 
served equivalent widths. The only possibility we have found 
is modifying the incident continuum such that = —2 
(F — 3). While inconsistent with the observed power law 
of r ~ 1 .75, this does yield equivalent widths consistent with 
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TABLE 1 

Comparison OF e quivale nt widths for PG 1211+143 (from 
I to a continuum driven wind. 



) MPARISON OF EQUIVALE N 

■Pounds et AL.I2003At ' 



Line 


Observed 
EW [eV]" 


Model EW [eV]'' 
Nu,o = 5 X 10^^ 


Model EW [eV]*"^ 
A^H = 5 X lO^^ 


Fe XXVI Lya 


95 ± 20 


7 X 10"'' 


31 


Fe XXV ls-3p 


45 ± 12 


2 X 10-5 


0.2 


Fe XXVI Ly/3 


45 ± 12 


1 X 10-3 


8 


S XVI Lya 


32 ± 12 


2 X 10-" 


1.4 


Mg XII Lya 


15 ±6 


6 X 10-^5 


0.4 



^EiTors listed correspond to 90% confidence intervals. 

^Equivalent widths are those found at s — 3ro, the infen'ed location of the sight-line through 
the outflow in PG 1211 + 143. 

'^The column for this calculation is held co nstant throughout the outflow, unlike all 
models, to compare with lPounds et alj|20033) . 



L/L-Edd > 1), completely stripping the outflowing 
gas, resulting in acceleration primarily by electron scat- 
tering. Continuum-driven winds can be launched by 
blackbody SEDs but such outflows do not have the ob- 
served high-ionization lines. 



3. Taking the particul ar example of the ou tflow observed 
in PG 1211H-143 (Pounds etal. 2003a), wind models 
display much smaller EWs than the observed EWs. 
Softening the X-ray power law to F ^ 3 can increase 
the EWs, however this requires both a non-standard 
AGN continuum and a constant-column disk wind. 
If pure continuum driving has difficulty reproducing the ob- 

previous . , .° , ■ , ■ 

servations, what process is accelerating the wind First, per- 



the observations for all lines except Fe XXV, but only for the 
wind with constant column of iVn = 5 x 10^^ cm-^ along 
the entire outflow. However, in the rather implausible case 
that the X-ray continuum is formed outside the wind (so that 
the wind sees a different X-ray continuum than is observed), 
such a model does reproduce the observed equivalent widths. 

As a final consideration, we note that for a large outflow 
with covering factor h = 0.8, P-Cygni profiles should be ob- 
served, but emission is not seen with the observed absorption 
lines. This may perhaps indicate a smaller covering fraction 
for the wind. 

4. SUMMARY 

This Letter addresses observational constraints on winds 
launched from near the central black hole of AGNs using pure 
continuum driving. This has become an increasingly signifi- 
cant topic very recently, with the observation of high-velocity, 
very highly ionized gas in several different QSOs, where the 
outflow is most likely too highly ionized to allow for signifi- 
cant line driving. As expected, continuum driving can accel- 
erate highly-ionized winds from near the central black hole to 
the observed velocities when L/L^dd > 1- But the models 
also indicate that: 

1. Changes in the initial density or shielding column do 
not affect the final velocity of the wind. Nor do changes 
in the illuminating continuum significantly change the 
final velocities in the wind, although moving to softer 
continua does significantly lower the ionization state. 

2. Continuum-driven winds, illuminated by an AGN SED 
and accelerating to the observed velocities quickly be- 
come very highly ionized (with ^ ~ 10^ — 10^ for 



haps the wind is not relativistic (see Kaspi 2004, who finds 
the low-energy X-ray spectrum of PG 1211H-143 to be con- 
sistent with a 3000 km/s outflow). It may also be possible 
that line-driving is playing a role, although such acceleration 
would have to occur only at the very base of the wind, where 
the ionization parameter is suitably low and where radiative 
launching would be difficult given the high column to the 
central source near the disk (but see Proaa & Kallman 2004). 
Another option for explaining these outflows is a multi-phase 
wind, although such a wind would increase the mass outflow 
rate beyond its already large value. 

Another possibility is that such high-velocity, highly- 
ionized winds could be launched hydromagnetically. The ad- 
dition of hydromagnetic driving would decouple the ioniza- 
tion state from the acceleration process; instead of requiring, 
as for continuum-driven flows, that the luminosity be iEdd 
or greater (which leads to over-ionization), the wind could 
possibly be accelerated by MHD forces that would allow a 
lower L/L-Edd and a lower ionization state. Indeed, a lower 
L/L-Edd may already be indicated by the larger A/bh (albeit 
with large uncertainty) found by Peterson et al. (2004). MHD 
wind models will be addressed in a subsequent paper, but we 
point out that the large mass outflow rates (> AfEdd) inferred 
in this source may be too high for an MHD outflow. Model 
predictions for winds driven by all of the above processes will 
be useful in determining which forces are capable of launch- 
ing these outflows and which model outflows are consistent 
with the observational signatures. 
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